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Abstract 

It is shown that the geometrical optics limit of the Maxwell equa- 
tions for certain nonlinear media with slow variation along one axis 
and particular dependence of dielectric constant on the frequency and 
fields gives rise to the dispersionless Veselov-Novikov equation for re- 
fractive index. It is demonstrated that the last one is amenable to the 
quasiclassical 9-dressing method. A connection is noted between ge- 
ometrical optics phenomena under consideration and quasiconformal 
mappings on the plane. 

PACS numbers: 02.30.Ik; 42.15.Dp 

Great variety of nonlinear phenomena from various fields of physics, ap- 
plied physics, mathematics and applied mathematics can be modeled by non- 
linear integrable equations PO-[H]- A subclass of such equations , the so-called 
dispersionless integrable equations, has attracted a particular interest during 
the last decade P]-[T7j. 

In this letter, we will show that the propagation of electromagnetic waves 
of high frequency in certain nonlinear media is governed by the dispersion- 
less Veselov-Novikov (dVN) equation. Namely, we will demonstrate that the 
Maxwell equations describing the propagation of waves in media character- 
ized by slow variation along the axis z and particular dependence of dielectric 
constant on frequency and fields in the limit of geometrical optics give rise 
to the dVN equations for refractive index. These equations provide us with 



integrable deformations of the plane eikonal equation which preserve, in par- 
ticular, the total "plane" squared refractive index J J n^dxdy. Under more 
special conditions one gets the dispersionless Kadomtsev-Petviashvili (dKP) 
equation. We will show that the dNV equations are treatable by the quasi- 
classical 9-dressing method developed recently in [HI- EDI- The quasiclassical 
5-dressing method reveals also a connection between the geometrical optics 
phenomena under consideration and the theory of quasiconformal mappings 
on the plane. 

We start with the Maxwell equations in the absence of sources (we put 
the velocity of the light c = 1) 

VAH-D = V-D = 

VAE + B = V-B = (1) 

and the material equations of a medium 

D = eE, B = fiU. (2) 

Here and below V = [dx, dy, dz), while ■ and A denote the scalar and vector 
product respectively. 

We will study a propagation of electromagnetic waves of the fixed, high 
frequency u, i.e. we will look for the solutions of the Maxwell equations of 
the form [21] 

E {x, y, z, t) = Eo {x, y, z) e-^'^*+^'^^(^'f'^) 

H {x, y, z, t) = Ho (x, y, z) (3) 

where Eq, Hq and the phase S (x, y, z) are certain functions. 
In addition we assume that the nonlinear medium is characterized 1) by the 
independence of the magnetic permeability //(x, y, z) on the frequency oj; 2) 
by the Cole-Cole dependence 
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eo + ^^^ Z2U^ 0<i/<- (4) 



of the dielectric constant e on u; (where e^^x^y^z) and tq are independent 
on uj and e is a function of coordinates and fields); 3) by slow variation of 
all quantities along the axis z such that ^ = uj^'^-^, where ^ is a "slow" 
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variable defined hy z = cu^^, i.e. 

S = S{x,y,^) 

Ho = Iioix,y,Lu''0- (5) 

Now, rewriting equations ((H) and Q as the second order differential equa- 
tions for the electric and magnetic fields, respectively, using Q and (0) and 
taking into account that at a; ^ oo 

^ = ^o + J, (6) 

where ei = (iro)'^'^ e, one gets in the leading u;^ order 

Sl + Sl = /xeo (7) 
while, in the next o;^"^'^ order, one obtains 

S^ = ifie^)K (8) 

Equation ((Tj) is the standard eikonal equation on the x, y plane with the 
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refractive index n{x,y,^) = (/ieo)^. Slow variation of S along 2;-axis is 
defined by equation (jH)), where fiei depends both on coordinates x,y,^ and 
fields. We also assume that the absorption effects are negligible, so /xe can 
be taken as real one . We would like to note that the first two features of the 
medium mentioned above are quite generic and they are valid for a variety of 
dielectric media P2|-|21I, while the particular type of slow variations along 
2;-axis given by (jH)), seems, did not attract attention before. 

Since the function ei is determined by the geometrical optics limit u oo 
and the time translation symmetry t ^ t + const of the Maxwell equations 
for the solutions of the type © is equivalent to the symmetry under trans- 
formations S ^ S + const, one concludes that the function ei for nonlinear 
media may depend on the coordinates x,y,^ and only on and Sy. Thus, 
the geometrical optics limit of the Maxwell equations in nonlinear media 
under consideration is governed by the equations 

Sl + S'y = n^x,y,0, (9) 
= 'fix,y,^;S.j:,Sy) (10) 
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where is some real-valued function. Introducing the variables z = x + iy^ 
z = X — iy (please, do not confuse it with the variable used in equation (jS))), 
one rewrites these equations as follows 

S,S, = u{z,z,0 (11) 
S^ = if {z,z;^, 3^,82) , (12) 

where u = An^. 

The compatibility of equations © and (fTUI) (or (fTT|). (fT^ ) imposes con- 
straints on the possible forms of the function (f, namely 

dip „ 

bz— + bz— + Uz'P +U2'P = u^, (13) 

where 

f' = -Q^{z,z;S^,Sz) , ip" = —{z,z;S,,Sz) . (14) 

Here we restrict ourself by functions which are polynomial in Sz, Sg and 
compatible with real-valuedness of S and u. For the simplest choice (p = 
ao{z, z, equation (fT3|) obviously gives ao = const, i.e. = 0, S = ao^ + 

S{z, z). For the linear function ip = aSz + (ySz + P + j3, one gets a = a{z), [3 = 
const, and 

= (au)^ + (au)- . (15) 

For the quadratic ip = aS^ + aS'^ + pSz + J^Sg + j + j, equations (fT^ and (fTTjl 
imply a = 0, P = P{z) and 7 = const, i.e. one gets the previous linear case. 

The cubic tp = aS'^ + aS^ + (3Sl + ]3Sl + -iS^ + ^S^ + 5 + 5 obeys 
equations (fT^ and (fTT|) if 

a = a{z), P = 0, 
7g = —azU — 3aUz, 

6 = const, (16) 

and one gets 

= + (7u)- . (17) 

In the particular case a = 1 and, consequently 7^ = —3uz, equation (fT7|) is 
nothing but the dispersionless Veselov-Novikov (dVN) equation introduced 

in [nii^ni- 

In a similar manner one can construct nonlinear equations which corre- 
spond to higher degree polynomials p:>. These higher degree cases apparently 
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become physically relevant for the phenomena with large values of and 
Sy. Thus, if we formally admit all possible degrees of and Sz in the right 
hand side of equation ()12j) . then one has an infinite family of nonlinear equa- 
tions, which may govern the ^-variations of the wave fronts and "refractive 
index" u. Since equation (fT^ should respects the symmetry S —S of 
the eikonal equation one readily concludes that only polynomials if of 
the form J2n=i '^nS^^'^ + J2n=i '^nSg^~^ , are admissible (the constant terms 
which have appeared in the cases (fT3j) . (fT7|) discussed above are, in fact, ir- 
relevant). In the case un = ^ one gets the dVN equation mentioned above 
{N = 2) and the so-called dVN hierarchy of nonlinear equations. The dVN 
equation has been introduced in jTH [20] as the dispersionless limit of the 
VN equation, which is the 2 + 1-dimensional integrable generalization of the 
famous Korteweg-de-Vries (KdV) equation PP-[3],[H]. 

Let us consider a more specific situation in which the propagation of elec- 
tromagnetic waves in the media discussed above exhibits also a slow variation 
along the axis y, namely dy = e9^, where 77 is a slow variable defined hj y = - 
and e is a small parameter. If one assumes that the functions S and u in the 
eikonal equation © have the following behavior at small e 

S = ^ - S {x,r],^) 

u(^x,^,C^ = ^ - q{x,ri,^) (18) 

and the polynomial (x, |, ^; S^, eS^^) has an appropriate behavior as e ^ 0, 
then in this limit, equations Q and (fTUI) are reduced to 

Sr, = Sl + q 

S^ = g:>(^x,'r],^]S^^ (19) 

where if is an odd order polynomial in S^. Equations (fT^ describe propa- 
gation of the quasi-plane wave fronts y = const in a medium with very large 
refractive index. Compatibility of equations (|19p gives rise to the well known 
dispersionless Kadomtsev-Petviashvili (dKP) equation = ^qq^ + ^d'^qrjr, 
and the whole dKP hierarchy. The dKP equation is rather well studied (see 
e.g. fn\ [T3] and reference therein). The KP equation itself represents the 
most known 2+1-dimensional integrable generalization of the KdV equation. 

The dVN and the dKP equations being relevant in particular situations 
of propagation of waves in certain nonlinear media have an advantage to 
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be solvable. Various methods have been applied to solve the dKP equation, 
including the quasiclassical 9-dressing method |IH|-|2D|- Here we will demon- 
strate that the dVN equation is amenable to this method too. 
The quasiclassical (9-dressing method is based on the nonlinear Beltrami 
equation 

S-, = W{X,X;Sx), (20) 

where A is the complex variable, S'a = |f and W is a certain function {d- 
data). To construct integrable equations one has to specify the domain G 
(in the complex plane C) of support for the function W {W = 0, A G C \ G) 
and look for solution of (pUj) in the form S = So + S, where the function 
Sq is analytic inside G, while S is analytic outside G [IE]- 1201 • To get the 
dVN equation and the whole dVN hierarchy, we choose G as the ring V = 
{AgC: i<|A|<a}, where a is an arbitrary real number (a > 1), and the 
function So in the form 

oo oo 

So = J2x^X'--' + Y,x^X~''^+\ (21) 

71=1 n=l 

We also assume that 

oo qioo) 

S = J:^; A-oo (22) 

n=0 

"normalization" and 

oo 
n=0 

Note that a function S (A, A, x„) has the properties ()2HI23|) if it obeys the 
constraints 

5(-A,-A) =-5(A,A), (24) 
S {X,X) = S (^l,\y (25) 

The constraint also implies that 

r.(oo) _ c>(0) 
'~'2n+l — '-'2n+l- 
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An important property of the nonlinear 9-problem (|2Up is that the derivatives 
of 5* with respect to any independent variable x„ obeys the linear Beltrami 
equation 




where W (A, A; (p) = ^ (A, A; 0). Equations has two basic properties, 
namely, 1) any differentiable function of a solution is again a solution; 2) 
under certain mild conditions on W, a solution which is equal to zero at 
certain point Aq G C, vanishes identically (Vekua's theorem) [23]. The use of 
these two properties allows us to construct algebraic equations of the form 

Denoting xi = z, choosing X2 = X2 = ^, taking into account that = 
A + Sz, Sz = A^^ + Sz, = (A'^ + ^) + S^, and using the properties of the 
linear Beltrami equation mentioned above, one gets 

SzSz = u, (27) 

= Si + si + VSz + VSz, (28) 

where 

u {z, z,0 = l + SiT\ V = -3S[f = -3d^\z. (29) 

Evaluating the terms of the order ^ in the both sides of equation (j28|) . one 
gets the dVN equation 

= -3 (tiSj^M^)^ - 3 {ud~\s)g ■ (30) 

Considering higher "times" 0:3, X4,... such that Xn = x„, one obtains the 
equations 

n 

Sx„ = Yl K^r-' + u^Sj"'-') , n = 3, 4, . . . (31) 

m=l 

Equations (jHH) together with (j77|) and provides us with the infinite dVN 
hierarchy of equations. 

Thus, the quasiclassical (9-dressing method allows us to treat the eikonal 
equation (j2ZI), equations and which describe its deformations and 
the corresponding dVN hierarchy for u, in a way similar to the dKP and 
d2DTL hierarchies PHJ-Eni- These integrable dVN deformations of the eikonal 
equation (j27p are quite special. In particular, they are characterized by the 
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existence of an infinite set of integral quantities (integrals of motion for dVN 
equation), which are preserved by these deformations. The simplest of them, 
for u such that m — >• as |z| — > cxd, is given by the total squared refractive 
index Ci = / J u {z, z,^) dz A dz (Ci^ = 0). 

Constraint (j25j) guarantees that the refractive index u is a real one. In- 
deed, taking the complex conjugation of equation (j77j) . using the differential 
consequences (with respect to z and z), of the above constraint and taking 
into account the independence of the l.h.s. of equation (j2Zj) on A, A, one gets 



U iXr 



(^n) — Sz (a, a, Xr^ Sz (a, a, Xn) — Sz ^ ^ ' ^ ' "^"^ ^ A ' A ' "'""^ 

(32) 

Moreover, this constraint implies that the function S is real-valued on the 
unit circle |A| = 1 {S (A, A) =S{\X), |A| = 1). 

Thus, the approach proposed provides us with the real- valued refractive index 
u {xn) and the phase function S together with their slow variations in the 
direction orthogonal to the plane y. Note that the refractive index u = 4n^ 
not necessary should be positive. For certain media the product /ueo can be 
negative as well PE] . 

The quasiclassical 9-dressing method provides also an effective method to 
construct exact solutions of the dVN equations as in the dKP case [THj-pU]. 
Exact solutions of the dVN equation and their applications in geometrical 
optics will be discussed in a separate paper. 

Finally we would like to note that solutions of the nonlinear Beltrami 
equation represent themself the so-called quasiconformal mappings (for 
dKP equation see ^3). In the our case we have quasiconformal map- 
pings of the ring V, given by the function S (A, A, x„) , with the boundary 
conditions and and constraints and (j2SI)- Thus the quasi- 
classical (9-dressing approach to the eikonal equation (|77|) and its deforma- 
tions by dVN hierarchy establishes a connection between the geometrical 
optics and the theory of the quasiconformal mappings on the plane (23 12H] ■ 
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